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ABSTRACT 

We combine new dust continuum observations of the edge-on spiral galaxy NGC 4565 in all 
Herschel/SPIRE (250, 350, 500 fim) wavebands, obtained as part of the Herschel Reference 
Survey, and a large set of ancillary data (Spitzer, SDSS, GALEX) to analyze its dust energy 
balance. We fit a radiative transfer model for the stars and dust to the optical maps with the 
fitting algorithm FitSKIRT. To account for the observed UV and mid-infrared emission, this 
initial model was supplemented with both obscured and unobscured star-forming regions. 
Even though these star-forming complexes provide an additional heating source for the dust, 
the far-infrared/submillimeter emission long wards of 100 /urn is underestimated by a factor of 
3-4. This inconsistency in the dust energy budget of NGC 4565 suggests that a sizable fraction 
(two-thirds) of the total dust reservoir (Md ~ 2.9 x 10 8 M ) consists of a clumpy distribution 
with no associated young stellar sources. The distribution of those dense dust clouds would be 
KJ{ in such a way that they remain unresolved in current far-infrared/submillimeter observations 

and hardly contribute to the attenuation at optical wavelengths. More than two-thirds of the 
dust heating in NGC 4565 is powered by the old stellar population, with localized embedded 
sources supplying the remaining dust heating in NGC 4565. The results from this detailed dust 
energy balance study in NGC 4565 is consistent with that of similar analyses of other edge-on 
spirals. 

Key words: radiative transfer - dust, extinction - galaxies: ISM - infrared: galaxies - galax- 
ies: NGC4565 
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1 INTRODUCTION 

Panchromatic radiative transfer 2D and 3D modelling of galaxies 
provide a powerful tool to analyze the characteristics of dust in 
galaxies (i.e. optical properties, distribution, dumpiness, etc.) in a 
self-consistent way (e.g. Xilouris et al. 1999; Popescu et al. 2000; 
Popescu et al. 2011; Alton et al. 2004; Bianchi 2008; Baes et al. 
2010; MacLachlan et al. 2011; De Looze et al. 2012; Holwerda et 
al. 2012). From optical observations, the properties and spatial dis- 



tribution of stars and dust can be constrained using a radiative trans- 
fer code to model the propagation of stellar light and its interaction 
with dust particles in a galaxy. In a second step, the dust emis- 
sion predicted from the radiative transfer simulations is compared 
to the observed thermal dust re-emission at infrared/submillimeter 
wavelengths. Such a complementary study imposes that dust fea- 
tures can easily be identified from optical as well as infrared ob- 
servations. This requirement has limited the number of galaxies for 
which detailed dust energy balance studies have been attempted in 
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the past. Edge-on spirals are considered ideal cases for those com- 
plementary studies since projection effects allow to resolve the dust 
distribution vertically (e.g. Misiriotis et al. 2001; Alton et al. 2004; 
Dasyra et al. 2005; Bianchi 2008; Baes et al. 2010; Popescu et al. 
2011; De Looze et al. 2012). Unfortunately, the large inclination 
angle impedes the characterization of substructures in the dust or- 
ganization and the identification of heavily obscured star-forming 
regions. Only the application of complex 3D models in conjunction 
with competitive fitting algorithms have been shown to success- 
fully characterize dumpiness and asymmetries in the stellar and 
dust distribution of highly inclined objects (e.g. Schechtman-Rook 
et al. 2012). 

Dust energy balance studies of individual edge-on galaxies re- 
veal an inconsistency between the predicted FIR/submm fluxes of 
radiative transfer models and the observed emission in those wave- 
bands (e.g. Popescu et al. 2000; Misiriotis et al. 2001; Alton et al. 
2004; Dasyra et al. 2005; Baes et al. 2010; Popescu et al. 2011; 
De Looze et al. 2012; Holwerda et al. 2012). Although radiative 
transfer models might successfully explain the observed optical at- 
tenuation, the modelled dust emission underestimates the observed 
thermal dust re-emission by a factor of 3-4. In order to reconcile the 
results of the radiative transfer models with the observations, two 
scenarios have been proposed: either a significant underestimation 
of the FIR/submm dust emissivity has been argued (e.g. Alton et al. 
2004; Dasyra et al. 2005; MacLachlan et al. 2011) or, alternatively, 
the distribution of a sizable fraction of dust in clumps or a second 
inner dust disk having a negligible attenuation on the bulk of the 
starlight (e.g. Popescu et al. 2000; Misiriotis et al. 2001; Bianchi 
2008; De Looze et al. 2012; MacLachlan et al. 2011). However, a 
general consensus on the origin of the energy balance problem has 
not yet been achieved. 

Until recently, observations of dust in galaxies were ham- 
pered by the poor resolution and limited wavelength coverage of 
infrared instrumentation. The first far-infrared space satellites (e.g. 
IRAS, ISO, Spitzer, Akari) only covered mid- and far-infrared 
wavelengths up to 200 /mi. Most ground-based submillimeter/mm 
instrumentation either only start operating from 850 /im onwards 
(e.g. LABOCA, IRAM) or suffer severely from the earth's opacity 
at shorter submm wavelengths (450 /xm, e.g. SABOCA, SCUBA). 
This coarse sampling of the Rayleigh- Jeans tail of dust emission in 
most galaxies induced poor constraints on the temperature, content 
and properties of dust in those objects. Recently, this large void in 
wavelength coverage was bridged thanks to the launch of the Her- 
schel Space Observatory (Pilbratt et al. 2010), offering the opportu- 
nity to trace the dust emission from galaxies in a broad wavelength 
domain ranging from 55 to 672 /xm. The high spatial resolution 
achieved by Herschel is providing a large collection of galaxies 
suitable for detailed dust energy balance analyses, supplemented 
with the ongoing development of powerful radiative transfer codes. 

At a distance of about 16.9 Mpc, the Needle Galaxy 
(NGC 4565) is one of the most nearby edge-on spiral galaxies cov- 
ering an angular scale of ^ 15' on the sky. Distance measurements 
for NGC 4565 range from 9 to 22 Mpc depending on the applied 
technique. We adopt a value of D ~ 16.9 Mpc throughout this pa- 
per, which is the average of the most recent distance measurements 
from / band surface brightness fluctuations as reported in Tonry et 
al. (2001) and Jensen et al. (2003). All values from the literature 
referred to in this paper were converted to the assumed distance of 
16.9 Mpc. NGC 4565 is a Sb spiral, classified as a Seyfert galaxy 
(Ho et al. 1997) and located in the Coma I cloud. Recently its clas- 
sification as an active galactic nucleus has been confirmed in the 
mid-infrared with the detection of the [Ne V] lines with Spitzer 



(Laine et al. 2010), which are among the strongest indicators of 
AGNs (see e.g. Tommasin et al. 2010). The stellar geometry in 
NGC 4565 has been studied extensively at optical and NIR wave- 
lengths (van der Kruit 1979; Hamabe et al. 1980; van der Kruit & 
Searle 1981; Jensen & Thuan 1982; Dettmar & Wielebinski 1986; 
Rice et al. 1996; Naeslund & Joersaeter 1997; Wu et al. 2002; Seth 
et al. 2005; Mosenkov et al. 2010; Kormendy & Barentine 2010). 
From studies of the absorption layer in NGC 4565, Ohta & Kodaira 

(1995) claim a dust distribution in a ring-like structure and an inner 
cut-off radius along the major axis at ~ 130". At radii smaller than 
the cut-off radius dust obscuration effects seem to be absent, ei- 
ther due to a deficiency of dust or the distribution of a central dust 
component preventing the identification of dust attenuation from 
optical data (Ohta & Kodaira 1995). Studies of residual J and K 
band images obtained after fitting the stellar component and sub- 
tracting those models from the near-infrared images in Rice et al. 

(1996) result in estimates for the dust scale length ~ 66" and scale 
height ~ 5" in NGC 4565, which corresponds to about 60 and 50% 
of the stellar disk in the K band. 

In addition to these optical/NIR studies of the dust component, 
the properties and geometry of the dust in NGC 4565 have been 
analyzed as well from infrared/submillimeter observations (Wain- 
scoat et al. 1987; Engargiola & Harper 1992; Rice et al. 1996; 
Neininger et al. 1996; Alton et al. 2004; Kormendy & Barentine 
2010; Laine et al. 2010). In the mid-infrared IRAC 8/xm waveband, 
Kormendy & Barentine (2010) and Laine et al. (2010) identify a 
ring of prominent poly cyclic aromatic hydrocarbon emission from 
dust at the same position as the molecular gas ring at radii of ~ 
80-100" (Sofue & Nakai 1994; Neininger et al. 1996). The first far- 
infrared study on dust emission in NGC 4565 based on IRAS data 
calculate a total infrared luminosity Lfir ~ 3 x 10 9 L (Wain- 
scoat et al. 1987), which compares well to the dust emission from 
our Galaxy. The dust emission in NGC 4565 however is more ex- 
tended compared to the dust distribution in our Galaxy (Wainscoat 
et al. 1987). Based on estimates of the infrared excess, Wainscoat 
et al. (1987) argue that only one quarter of the dust reservoir in 
NGC 4565 is heated by embedded stars in molecular clouds with 
the remainder of the dust heating power provided by the diffuse 
interstellar radiation field. 

Observations at FIR wavelengths (100, 160, 200 jim) with the 
Kuiper Airborne Observatory (Engargiola & Harper 1992) find a 
two-component dust model with a cold dust reservoir (Td ~ 20 K) 
distributed in an exponential disk and a warmer dust component 
(Td ~ 30 K) spatially coinciding with a bisymmetric spiral pattern 
in NGC 4565. In spatial correlation with the more extended ring of 
Hi gas, Neininger et al. (1996) also found a colder dust component 
(~ 15 K) in the outskirts of the galaxy. A second dust reservoir (Td 
~ 18 K) was found in the center of NGC 4565 and in the molecular 
gas ring at radii of ~ 80-100 ". The central concentration of both 
dust and molecular gas (Neininger et al. 1996) contradicts the sug- 
gested absence of dust in the centre of the galaxy (Ohta & Kodaira 
1995). 

Alton et al. (2004) constrain the amount of dust in NGC 4565 
from a complementary study of the attenuation properties of dust 
in the optical and the thermal dust re-emission at 1.2 mm. From 
optical constraints, they construct a radiative transfer model for the 
stars and dust in NGC 4565. Upon comparison of the visual opti- 
cal depth in their radiative transfer model with thermal continuum 
radiation in the 1.2 mm waveband, Alton et al. (2004) infer a dust 
emissivity at 1.2 mm which is 1.5 times higher than the bench- 
mark, semi-empirical model from Draine & Lee (1984) used in 
the radiative transfer simulation. From their combined study on the 
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dust emissivity in NGC4565 and two other nearby spirals, they 
argue that coagulation of well-ordered dust crystalline grains into 
amorphous particles in high-density environments can influence the 
emissivity of dust at submm/mm wavelengths. 

The first submillimeter observations of NGC 4565 were ob- 
tained during the first test flight of the 2 m Balloon-borne Large 
Aperture Submillimeter Telescope (BLAST) in 2005 (Wiebe et 
al. 2009), covering the galaxy in the same wavebands (250, 350, 
500 /xm) as offered by the SPIRE instrument onboard the Herschel 
Space Observatory. Although the resolution in their final images 
was degraded by the large point- spread functions (PSF) (~ 186- 
189", Truch et al. 2008), Wiebe et al. (2009) could estimate the 
scale length of the dust lane Kr ~ 150" from the beam-convolved 
image at 250 /xm, which is about twice the value obtained from 
optical studies (Rice et al. 1996). A modified black-body fit to 
the BLAST fluxes and additional constraints from ISO (170 /mi, 
Stickel et al. 2004) and IRAS (12, 25, 60, 100 /mi, Sanders et al. 
2003) observations predict a dust reservoir of mass Md ^3 x 10 8 
Mq (scaled to a distance of 16.9 Mpc) at a temperature Td ~ 16 K 
for an emissivity index f3 = 2 and a Draine & Li (2007) dust grain 
composition (Wiebe et al. 2009). 

NGC 4565 was recently observed in all SPIRE wavebands 
(250, 350 and 500 /xm) as part of the Herschel Reference Survey 
(HRS, Boselli et al. 2010). In this work, we will analyze the dust 
characteristics from a spatially resolved study of the Herschel in- 
frared/submillimeter photometry and a wealth of ancillary observa- 
tions. Section 2 reviews the observing and data reduction strategy 
for the Herschel observations and the set of ancillary data used in 
this study. Section 3 presents a detailed study of the dust energy 
balance in NGC 4565 with a brief introduction to the applied radia- 
tive transfer code SKIRT and fitting algorithm FitSKIRT (Section 
3.1) and a gradual construction of a self-consistent model that can 
account for the observed quantities for NGC 4565 across the multi- 
wavelength spectrum (Section 3.2). In Section 4, the results of our 
radiative transfer modelling procedure are discussed and compared 
to the dust energy balance, gas-to-dust ratio and dust heating mech- 
anisms in other nearby galaxies. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Herschel data (250-500 micron) 

NGC 4565 was observed on the 31st of December 2009 (ObsID 
0x50002cd8) as part of the Herschel Reference Survey (HRS, 
Boselli et al. 2010), a Herschel Guaranteed Time Key Program 
observing 323 galaxies in the nearby Universe spanning a wide 
range in morphological type and environment. Dust continuum 
maps were obtained in all SPIRE wavebands (centered at 250, 350 
and 500 /mi), covering the galaxy in three cross-linked scans (nom- 
inal and orthogonal) at medium scan speed (30"/s). The final maps 
cover a square area of 30' x 30', corresponding to an extent of at 
least 1.5 times the optical isophotal diameter D 2 5 in NGC 4565. 

The SPIRE data were processed up to Level- 1 in the Herschel 
Interactive Processing Environment (HIPE, version 4.0.1367, Ott 
2010) with the standard script adapted from the official pipeline 
(POF5 pipeline.py, dated 8 Jun 2010) as provided by the SPIRE In- 
strument Control Centre (ICC). Instead of the ICC default settings 
for deglitching, we used the waveletDeglitcher, which was adjusted 
to mask the sample following a glitch. The waveletDeglitcher algo- 
rithm was applied a second time after the flux calibration, since 
this was shown to significantly improve the removal of remaining 



Table 1. Panchromatic overview of flux densities used in the dust energy 
balance study of NGC 4565 



Filter 


A (/mi) 


?v (Jy) 


Kei 


FUV 


0.15 


0.008 ± 0.001 


1 


NUV 


0.23 


0.015 ± 0.002 


1 


SDSSu 


0.36 


0.07 ± 0.01 


2 


SDSSg 


0.47 


0.33 ± 0.03 


1 


SDSSr 


0.62 


0.70 ± 0.07 


1 


SDSSi 


0.75 


1.08 ±0.10 


1 


SDSSz 


0.89 


1.45 ±0.14 


2 


2MASS J 


1.25 


2.18 ±0.03 


3 


2MASS H 


1.65 


3.05 ±2.51 


3 


2MASS K 


2.20 


2.51 ±0.04 


3 


WISE 1 


3.4 


1.54 ±0.04 


4 


IRAC 3.6 


3.6 


1.39 ±0.14 


5 


IRAC 4.5 


4.5 


0.91 ±0.09 


5 


WISE 2 


4.6 


0.90 ± 0.03 


4 


IRAC 5.8 


5.8 


1.17 ±0.14 


5 


IRAC 8.0 


8.0 


1.85 ±0.19 


5 


WISE 3 


12.1 


2.20 ± 0.05 


4 


IRAS 12 


12 


1.53 ±0.23 


6 


WISE 4 


22.2 


1.60 ± 0.16 


4 


MIPS 24 


24 


1.65 ±0.07 


6 


IRAS 25 


25 


1.70 ±0.26 


6 


IRAS 60 


60 


9.83 ± 1.47 


6 


MIPS 70 


70 


19.3 ± 1.9 


6 


IRAS 100 


100 


47.23 ± 7.08 


6 


KUIPER 100 


100 


52.00 ± 10.40 


8 


MIPS 160 


160 


86 ± 10 


7 


KUIPER 160 


160 


62.00 ± 12.40 


8 


KUIPER 200 


200 


55.00 ± 11.00 


8 


SPIRE 250 


250 


63.24 ± 1.18 


9 


SPIRE 350 


350 


31.35 ±0.69 


9 


SPIRE 500 


500 


13.14 ±0.29 


9 


IRAM 1200 


1200 


1.00 ±0.10 


10 



References: (1) Cortese et al. 2012b; (2) Cortese (priv. comm.); (3) Jarrett 
et al. 2003; (4) this paper, see Section 2.3; (5) Ciesla et al. (in prep.); (6) 
Rice et al. 1988; (7) Bendo et al. 2012; (8) Engargiola & Harper 1992; (9) 
Ciesla et al. 2012; (10) Neininger et al. 1996 

glitches for the version of the code at that time (currently it is only 
run once). Instead of the pipeline default temperature drift correc- 
tion and median baseline subtraction, we applied a custom method 
(BriGAdE; Smith et al. in prep.) to remove the temperature drift 
and bring all bolometers to the same level. Final SPIRE maps (see 
Figure 1 , last three bottom panels) were created with the naive map- 
per in HIPE with pixel sizes of 6", 8", and 12" at 250, 350, and 
500 /mi, respectively. The FWHM of the SPIRE beams have sizes 
of 18.2", 24.5", and 36.0" 1 at 250, 350, and 500 /mi (for pixel 
sizes of 6, 8 and 12"), respectively. To update our flux densities to 
the latest v8 calibration product, the 350 \m\ images are multiplied 
by a factor of 1.0067 (SPIRE Observers' Manual 2011). Since the 
calibration in the standard data reduction pipeline was optimized 
for point sources, correction factors (0.9828, 0.9834 and 0.9710 at 
250, 350 and 500 /xm) are applied to convert the K4 colour cor- 
rection factors from point source to extended source calibration 
(SPIRE Observers' Manual 2011). Appropriate colour correction 
factors were determined from a single component modified black- 
body fit to the Rayleigh- Jeans tail of dust emission (constrained 

1 Retrieved from document available at http://herschel.esac.esa.int/ 
twiki/pub/Public/SpireCalibrationWeb/beam_release_note_v 1 - 1 .pdf 
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by fluxes from MIPS 160 /im and all SPIRE wavebands) for values 
f3 = 1.0, 1.5 and 2.0 (i.e. the range of spectral indices derived for the 
HRS sample galaxies, Boselli et al. 2012). With the best fit obtained 
for values ft = 2 and Td = 17 K, we apply colour correction factors 
(1.019,1.009,1.021) to the 250, 350 and 500 /mi images (SPIRE 
Observers' Manual 2011), respectively, to correct for the shape of 
the response function in every filter. An uncertainty in the calibra- 
tion of 7% was assumed (Swinyard et al. 2010; SPIRE Observers' 
Manual 2011). 

From aperture photometry within an elliptical aperture (20 ; x 
2.7') flux densities of 62.06, 31.07 and 12.87 Jy for NGC 4565 were 
obtained at 250, 350 and 500 /mi, respectively (Ciesla et al. 2012). 
Since those fluxes were corrected for extended source calibration, 
we only need to correct those values with appropriate colour cor- 
rections factors (1.019, 1.009, 1.021). This results in final flux den- 
sities of 63.24, 31.35 and 13.14 Jy at 250, 350 and 500 /im, re- 
spectively. The large ancillary dataset of far-infrared/submillimeter 
observations allows us to compare those SPIRE fluxes to flux mea- 
surements from other space or airborne missions in overlapping 
wavebands. Wiebe et al. (2009) report fluxes obtained with BLAST 
at 250 /mi (37.2 ± 4.5 Jy), 350 /mi (21.0 ± 2.1 Jy ) and 500 /mi 
(9.8 ± 0.9 Jy), which are 40, 32 and 24 % lower with respect to 
SPIRE flux measurements (see Table 1). BLAST observations of 
NGC 4565 covered an area of ~ 0.4 deg 2 centered on the galaxy, 
sufficient to detect the dust emission from NGC 4565 in those 
wavebands. Considering that BLAST fluxes for NGC 4565 were 
measured from data taken during the first BLAST test flight in 
2005, the images were degraded by the large point- spread function 
of BLAST05 (FWHM of - 186-189", Truch et al. 2008). Prior to 
aperture photometry, the images were deconvolved following the 
method outlined in Chapin et al. (2008). This beam deconvolution 
might be the cause of the lower BLAST fluxes compared to SPIRE 
measurements. Alternatively, calibration issues might be responsi- 
ble for the small off-set between BLAST and SPIRE. Also Planck 
measurements (Planck Collaboration 2011) at 350 (24.35 ± 1.43 
Jy) and 550 /xm (8.12 ± 0.48 Jy) are about 20% below our SPIRE 
measurements. The lower Planck values can easily be explained 
from the radius applied for aperture photometry on point sources in 
their first release catalog (Planck Collaboration 2011). Indeed, an 
aperture with radius fixed to the FWHM of Planck at 350 /im (~ 
4.23') does not cover the entire galaxy's emission in that waveband 
(see the upper right panel of Figure 10 in Ciesla et al. 2012). Per- 
forming aperture photometry within the Planck radius at 350 /im 
on the SPIRE images results in a flux density S u ^ 25.5 ± 0.9 Jy 
(Ciesla et al. 2012), which is in close agreement with the reported 
Planck value (Planck Collaboration 201 1). 

Figure 1 (last three bottom panels) shows the Herschel maps 
for NGC 4565 in the SPIRE 250, 350 and 500 /xm wavebands. 
Weak warping signatures at the edges of the dust disk can be 
identified in all Herschel dust maps, with the warp on the north- 
western side of the galaxy's disk being more pronounced. These 
warp features were earlier identified in NGC 4565 from optical data 
(e.g. van der Kruit 1979; Jensen & Thuan 1982; Naeslund & Joer- 
saeter 1997), Hi (Sancisi 1976; Ruben 1991) and dust continuum 
(Neininger et al. 1996) observations and are most likely the relicts 
of a tidal interaction with its neighboring galaxy NGC 4562 (Jensen 
& Thuan 1982). 

Aside from the warps on the edges of the galaxy, three peaks in 
the dust emission can be identified in the SPIRE 250 and 350 /mi 
images (see Figure 1, second and third bottom panels). The two 
brightest emission peaks originate from a narrow ring structure cen- 
tered at a galacto-centric radius of 80-100", coinciding with a ring 



of molecular gas in NGC 4565 (Sofue & Nakai 1994; Neininger 
et al. 1996). In the center of NGC 4565 the emission peak corre- 
sponds to a central concentration of dust, which was already iden- 
tified from 1.2 mm observations (Neininger et al. 1996). Due to the 
growing size of the SPIRE beam towards longer wavelengths, the 
three peaks are more difficult to distinguish in the SPIRE 500 /xm 
image (see Figure 1, fourth bottom panel), only showing a small 
increase in emission at the edges of the narrow dust ring. 

Figure 2 compares the distribution of dust in the SPIRE 
250 /im image with the gas distribution from integrated Hi and CO 
line intensities (Neininger et al. 1996), the total gas column den- 
sity and dust emission profiles along the major-axis from 24 /xm 
and 1.2 mm dust continuum observations, smoothed to a common 
resolution of 20". Based on these major axis profiles, we can iden- 
tify a depression up to the 80" radius from the center in the Hi 
profile and a more or less constant column density profile in the 
outskirts of the galaxy from that radius onwards. The CO profile in 
NGC 4565 indicates the presence of a large molecular gas reservoir 
in the inner regions, distributed in a central disk and narrow ring 
with peak density around 80", where also the inner depression in 
Hi density ends. The SPIRE 250 /mi dust emission profile along 
the major axis correlates well with the molecular gas density pro- 
file (see also Wainscoat et al. 1987, Rice et al. 1996, Kormendy & 
Barentine 2010 and Laine et al. 2010) as well as the column den- 
sity of the combined Hi and H2 gaseous components. The latter 
column density profile seems dominated by the peaks in the molec- 
ular gas component, rather than the smooth Hi distribution across 
the plane of the galaxy. In contrast to the CO profile reported in 
Neininger et al. (1996), the CO observations from Sofue & Nakai 
(1994) identified a more asymmetric CO distribution with the ma- 
jority of the molecular gas reservoir residing on the northwest side 
of the galaxy. The dust emission from the colder dust reservoir in 
the SPIRE 250 /xm waveband however does not show a tendency 
for asymmetry. 

2.2 Spitzer data (3.6-160 micron) 

NGC 4565 was observed with the IRAC and MIPS instruments on- 
board the Spitzer Space Observatory as part of the programs Brown 
Dwarf Galaxy Haloes (PI: Giovanni Fazio; AOR key: 3627776) 
and The formation of dust lanes in nearby, edge-on disk galaxies 
(PI: Roelof S. de Jong; AOR key 14481408), respectively. IRAC 
(3.6, 4.5, 5.8, 8.0 /mi) and MIPS (24, 70, 160 /mi) data were re- 
trieved from the Spitzer archive and reduced following the proce- 
dure outlined in Ciesla et al.(in prep.) and Bendo et al. (2012). Al- 
though Spitzer IRAC data of NGC 4565 were already used in the 
analyses of Kormendy & Barentine (2010) and Laine et al. (2010), 
the maps presented in this work were processed independently. Fi- 
nal IRAC images have pixel sizes of 0.6", whereas MIPS images 
were processed to final maps with pixel sizes of 1.5", 4.5" and 9.0" 
at 24, 70 and 160 /mi, respectively. The IRAC beam has an almost 
uniform size in all wavebands (mean FWHM ~ 1.7", 1.7", 1.9" 
and 2.0" at 3.6, 4.5, 5.8 and 8 /mi, Fazio et al. 2004). The FWHM 
of the MIPS beam varies from 6" at 24 /mi, 18" at 70 /im to 38" at 
160 /mi (Engelbracht et al. 2007; Gordon et al. 2007; Stansberry et 
al. 2007). 

To determine flux densities for NGC 4565 from the processed 
IRAC images, stars in the surrounding field were masked to ob- 
tain realistic estimates for the background emission (Ciesla et al. in 
prep.). Once the images were background subtracted, flux measure- 
ments were derived within elliptical aperture matching the shape of 
the galaxy and encompassing all emission (see Table 1 and Ciesla 



The dust energy balance in the edge-on spiral galaxy NGC4565 5 



0.000 0.025 0.05 
MJy sr 1 



0.000 0.025 0.05 
MJy sr" 1 



0.0 0.5 1.0 

MJy sr" 1 



0.0 0.7 1.4 
MJy sr" 1 



0.0 0.7 1.4 



1.0 2.0 



0.0 1.5 3.0 



/ 1 / I / 



5 10 

MJy sr" 1 



IRAC 3.6 micron 



5 10 

MJy sr" 1 



10 5 

MJy sr" 1 



5 10 

MJy sr" 1 



45 90 

MJy sr" 1 



IRAC 5.8 micron 



MIPS 24 micron 



MIPS 70 micron 



65 130 
MJy sr" 1 




130 260 
MJy sr" 1 




MIPS 160 micron 



12:36.35 
Right Ascension (J2000) 



SPIRE 250 micron 



12:36.35 
Right Ascension (J2000) 



SPIRE 350 micron 



12:36.35 
Right Ascension (J2000) 



SPIRE 500 micron 



12:36.35 
Right Ascension (J2000) 



Figure 1. The stellar and dust emission from NGC4565 as observed from GALEX FUV and NUV, SDSS ugriz, Spitzer IRAC 3.6, 5.8 and 8.0 /mi, 
Spitzer MIPS 24, 70 and 160/im and Herschel SPIRE 250, 350 and 500 /im images. The IRAC 4.5 fj,m map is not explicitly shown here because of the great 
resemblance with the IRAC 3.6 fim image. The beam sizes of all instruments in the infrared/submillimeter wavebands have been indicated as white circles in 
the bottom left corner of the corresponding panels. 



et al. in prep.). Similarly, the identification and removal of fore- and 
background emission from MIPS data required some efforts before 
fluxes could be determined from the background subtracted MIPS 
images (Bendo et al. 2012, see Table 1). Uncertainties for IRAC 
and MIPS photometry measurements include uncertainties owing 
to the calibration, background noise and map making. For IRAC, 
the calibration is assumed to be accurate within 1.8, 1.9, 2.0 and 
2.1 % in IRAC 3.6, 4.8, 5.8 and 8.0 /xm wavebands (Reach et al. 
2005). For MIPS, calibration uncertainties of 4, 10 and 12 % were 
reported for the 24, 70 and 160 /xm wavebands (Engelbracht et al. 
2007; Gordon et al. 2007; Stansberry et al. 2007). 

MIPS fluxes at 24 and 70 /xm compare reasonably well with 



IRAS observations of NGC4565 at 25 and 60 /xm (see Table 1). 
The longer wavelength MIPS observations at 160 /xm are further- 
more consistent within 30% with data from the Kuiper Airborne 
Observatory (KAO) in the same waveband (see Table 1). 

Figure 1 (last panel on the second row) shows the stellar emis- 
sion in the IRAC 3.6 /xm image. Dust emission from IRAC (5.8 and 
8 /xm) and MIPS observations are presented in Figure 1 (third row 
and first panel on fourth row). The emission in the IRAC 3.6 and 
4.5 /xm wavebands is dominated by the stellar emission from the 
bulge in NGC4565. The stellar emission from the bulge is also 
discernible in the IRAC 5.8 /xm image. At 5.8 and 8 /xm, emis- 
sion from polycyclic aromatic hydrocarbons (PAHs) and very small 
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Figure 2. Major-axis profiles of the dust emission in MIPS 24 fim (red solid curve), SPIRE 250 /im (cyan dotted curve), IRAM 1.2 mm (yellow dashed 
curve) wavebands, the gas distribution from the Hi (black dashed-dotted curve) and CO (green solid curve) temperature brightnesses and the column density 
of gas combining the atomic and molecular components in NGC4565 (blue dashed curve). The Hi column density could be obtained directly from the 
Hi temperature brightness (N^i = 1.822 x 10 18 cm -2 jT^iy)^). For the calculation of the H2 column density from the CO temperature brightness 
(Nh 2 = Xco fTb(v)dv), we assumed a factor Xqo representative for the H2-to-CO fraction in our own Galaxy {Xqo ~ 2.3 x 10 20 cm -2 (K km 
s -1 ) -1 , Strong et al. 1988). All data were smoothed to a resolution of 20 /x . For ease of reference, all profiles were normalized by their maximum value. Those 
normalized profiles were then given an offset in the vertical direction to facilitate a mutual comparison. 



grains (VSGs) distributed in a ring centered at 80-1 00" dominates. 
The presence of a hot dust component heated predominantly by 
star formation is also clearly detected in the MIPS 24 /im image. 
Besides this ring of hot dust, the central region of NGC4565 is 
also a prominent emission source in the IRAC 5.8, 8 /im and MIPS 
24 /im wavebands. This central dust concentration in NGC 4565 
is most likely heated by the AGN, since the inner 80" region of 
NGC 4565 seems to lack any star formation activity (see Section 
2.4). Owing to the coarser resolution at 70 and 160/im, the dust 
emission features are not as easily identified compared to shorter 
wavelength data. To better contrast the structures in the far-infrared 
dust emission, we overlaid the contours of dust emission on those 
far-infrared MIPS data. The contours show three peaks in the far- 
infrared dust emission of NGC 4565, similar to the peaks identified 
from the Herschel data, with prominent emission from a ring at 
80-100" and a central dust concentration. Also the dust warps are 
manifested on the edges of the dust disk, in particular on the north- 
western side, in both IRAC (5.8, 8 /im) and all MIPS wavebands. 

In a similar way as for the SPIRE observations at 250 /im, the 
resolution in the MIPS 24 /im map was degraded to 20", allow- 
ing an immediate comparison with the distribution of atomic and 
molecular gas along the major axis of NGC 4565 (see Figure 2). 
With the shape of the PSF in the MIPS 24 /im filter deviating from 
a Gaussian distribution, appropriate kernels were constructed fol- 
lowing the procedure outlined in Bendo et al. (2012) to degrade 
the resolution in the MIPS 24 /im image. With the 24 /im emis- 
sion tracing the warm dust component in NGC 4565, we find a spa- 
tial correspondence with the dust heated by young stellar objects 
in NGC 4565 and the molecular gas reservoir, which provides the 
birth material for this star formation. The resemblance between the 
MIPS 24 /im, SPIRE 250 /im and IRAM 1.2 mm dust profiles sug- 
gests that the cold and warm dust component are well mixed in 
the interstellar medium of NGC 4565. However, the asymmetry in 
CO observations from Sofue & Nakai (1994) does become evident 
in the dust emission profiles of the warmer dust component. This 
asymmetry is most likely the consequence of spiral arm structures 



in the plane of the galaxy. Similar asymmetric absorption features 
could also be identified in optical images (Hamabe et al. 1980; Rice 
et al. 1996) and were assigned to a spiral arm structure at position 
angle -135°. 

2.3 Other infrared/(sub)millimeter data 

To refine the sampling of the spectral energy distribution, we in- 
clude photometric data, other than those obtained by Herschel and 
Spitzer, ranging from the infrared to millimeter wavelength do- 
main. At near- and mid-infrared wavelengths, we complement our 
dataset of observations with flux measurements from the Wide-field 
Infrared Survey Explorer (WISE, Wright et al. 2010). Images be- 
longing to the WISE all-sky survey at 3.4, 4.6, 12.1 and 22.2 /im 
were retrieved from the NASA/IPAC Infrared Science Archive. To 
measure fluxes from those images, all stars and background sources 
were masked in the field of the galaxy. Since the galaxy is located 
at the edge of the map, a determination of the background level 
from an annulus around the galaxy was impossible. A background 
value was determined from averaging the flux over 100 apertures 
taken randomly within a 30' field in this masked image. The size of 
the apertures was adapted to the FWHM of the PSF in every wave- 
band with a diameter of 4 x FWHM. The FHWM in the differ- 
ent WISE bands varies from 6.1", 6.4", 6.5" at 3.4, 4.6, 12.1 /im 
up to 12" in the longest 22.2 /im waveband (Wright et al. 2010). 
The mean background value was subtracted from the masked im- 
age. From this masked and background- subtracted image, we de- 
termined fluxes in all wavebands from aperture photometry. With 
the stellar bulge being a prominent emission feature at near-infrared 
wavelengths and the mid-infrared images having a dominant contri- 
bution from PAHs and very small grains in the star-forming disk of 
NGC 4565, the size of the aperture was adjusted to fit the size of the 
galaxy emitting in each waveband. Fluxes of the stellar emission 
from WISE at 3.4 and 4.6 /im are consistent with IRAC flux mea- 
surements at 3.6 and 4.5 /im (see Table 1). Also the WISE emis- 
sion in the longer wavelength bands corresponds well with mea- 



The dust energy balance in the edge-on spiral galaxy NGC4565 7 



surements by IRAS and MIPS at similar wavelengths (see Table 
1). 

IRAS flux measurements at 12, 25, 60 and 100 /xm were 
adopted from the Catalog of IRAS observations of large optical 
galaxies (Rice et al. 1988). Also the Kuiper Space Airborne Obser- 
vatory observed NGC4565 (Engargiola & Harper 1992) in wave- 
bands centered at 100, 160 and 200 /im (see Table 1). ISO observa- 
tions at 170 /xm (32.5 ± 4.9 Jy, Stickel et al. 2004) were not con- 
sidered here due to aperture effects (see also Wiebe et al. 2009). 
With Akari measurements at 65, 90, 140 and 160 /xm (1.17 ± 0.12 
Jy, 3.99 ± 0.22 Jy, 11.85 ± 1.07 Jy and 12.58 ± 2.39 Jy, respec- 
tively) being up to one order of magnitude lower in comparison 
to results from other airborne and space telescope facilities (IRAS, 
ISO, MIPS, KAO) in overlapping wavebands, we decided not to 
take into account Akari flux measurements for the analysis in this 
paper. 

At millimeter wavelengths, IRAM observations at 1.2 mm 
constrain the emission from cold dust in those wavebands. Alton et 
al. (2004) reported a flux density of 55 mJy when integrating along 
the major axis profile of the 1.2 mm emission for a 20" beam. As- 
suming the 1 .2 mm emission is spatially distributed similar to that 
at 250 /xm, we estimate a total flux density of 99 mJy at 1.2 mm for 
NGC4565 (see Table 1). 

2.4 GALEX and optical data (0.15-0.89 micron) 

GALEX FUVINUV and optical SDSS ugriz data for NGC 4565 
were recovered from the GALEX and SDSS archives and reduced 
following the instructions outlined in Cortese et al. (2012b). Fi- 
nal maps were reduced to obtain pixel sizes of 1.5" and 0.4" for 
GALEX and SDSS observations, respectively. GALEX FUV and 
NUV channels are characterized by a FHWM of 4.2" and 5.3", re- 
spectively (Morrissey et al. 2007), whereas the shape of the PSF in 
SDSS maps depends on the seeing during sky exposures. GALEX 
images for NGC 4565 were obtained in a single run of 1693 s expo- 
sure time, with SDSS images covered in 54 s of observing time. The 
first two panels on the top row of Figure 1 show the GALEX FUV 
and NUV images, respectively. Both GALEX images clearly show 
the low level of star formation activity in the inner most regions in 
NGC 4565. This lack of direct UV emission from young stars in 
the inner 80" region in NGC 4565 was already identified by En- 
gargiola & Harper (1992) and Ohta & Kodaira (1995), relating this 
deficiency in star formation to a bar structure in the central regions. 

Figure 1 presents the SDSS u, g, r, i and z band images in 
the last two top panels and first three panels on the second row. 
The main emission features in the optical SDSS images are the 
stellar disk and bulge. At shorter optical wavelengths (in partic- 
ular u band), a large fraction of the stellar light is extinguished 
along the line-of- sight by the prominent dust lane. Towards longer 
wavelengths, the dust obscuration diminishes and also the boxy 
shape of the bulge becomes more evident. This boxy shaped bulge 
is thought to result from a bar in the inner region of NGC 4565 
viewed edge-on (Combes & Sanders 1981) and is therefore con- 
sidered to be a pseudo-bulge. A genuine bulge component is also 
present in NGC 4565, but its spatial extent and intensity is infe- 
rior to the dominant emission from the pseudo-bulge (Kormendy & 
Barentine 2010). 

Both from GALEX and SDSS images, the warping structures 
on the southeastern and northwestern sides of the galaxy's disk can 
be discerned. Especially at UV wavelengths, the warps are promi- 
nent emission features, indicative for the ongoing star formation 
in these tidally disrupted parts of the interstellar medium (ISM) 



in NGC 4565. Flux densities for the stellar emission in UV/optical 
wavebands were adopted from Cortese et al. (2012b) for GALEX 
FUVINUV and SDSS gri bands (see Table 1). Fluxes for the 
remaining SDSS bands (u, z) were obtained in an homogeneous 
way from aperture photometry similar to the procedure for the other 
SDSS bands. 



3 DUST ENERGY BALANCE IN NGC 4565 

3.1 SKIRT and FitSKIRT 

SKIRT (Baes et al. 2003, 201 1) is a 3D Monte Carlo radiative trans- 
fer code designed to model the absorption, scattering and thermal 
re-emission of dust in a variety of environments: circumstellar disks 
(Vidal & Baes 2007; Vidal et al. 2011), clumpy tori around active 
galactic nuclei (Stalevski et al. 2012) and different galaxy types 
(Baes et al. 2010; De Looze et al. 2010; Gadotti et al. 2010; Gomez 
et al. 2010). 

FitSKIRT (De Geyter et al. 2012) is a fitting routine that com- 
bines the output of SKIRT with a genetic algorithm optimization 
library to obtain a best fitting model for the stellar and dust com- 
ponents in a galaxy. In this manner, the best fitting parameters for 
a variety of models and stellar and dust geometries can be derived 
when providing an input image from which the stellar emission and 
dust obscuration effects can be modelled. 

3.2 Radiative transfer modeling 

3.2.1 Model 1 

A first model accounting for the observed optical properties of stars 
and dust in NGC 4565 is obtained from the FitSKIRT fitting algo- 
rithm. FitSKIRT was applied to determine the best fitting model 
parameters using our SDSS g band image (see Figure 1, last top 
panel) as a reference image. We choose the g-band image for the 
fitting procedure, since the emission in this waveband is dominated 
by the old stellar population. At shorter wavelengths, the contribu- 
tion from young stars increases rapidly. The g band image is also 
optimal to perceive the attenuation effects of dust, which dimin- 
ishes at longer wavelengths. Prior to the FitSKIRT fitting algorithm, 
the stellar warps on both ends of the galactic disk were masked in 
this image. 

The input model for the stellar component consists of three 
components. The first component is an exponential disk with scale 
length Kr and height h z to fit the disk component 
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po exp 
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The second component is a flattened Sersic geometry for the outer 
bulge component. It is modelled as 
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where m = y R? + ^ is the spheroidal radius (q is the intrinsic 
flattening or axial ratio) and <S n (s) represents the dimensionless 
3D spatial density of a model that deprojects to a projected surface 
brightness profile 
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where 6„ is a dimensionless constant that can be approximated by 



8 /. De Looze et al. 



k + + o KK 1 6 K 2 + haIiIk s » as derived 

3 405n 25515n^ 1148175n d ' 

by Ciotti & Bertin (1999). Finally, the third component consists of 
an inner bulge described as a Sersic function (same as Equation 2 
but with m = VR 2 + z 2 ) with a Sersic index n ^ 1.4 and a scale 
length r e ff ~ 1.3", similar to the values reported in Kormendy & 
Barentine (2010). The luminosity of the inner bulge was scaled to 
a fixed value Ly ~ 3.6 x 10 4 Lv,©, optimized to reproduce the g 
band minor- axis profile (see Figure 3, last bottom panel). 

For the spectral energy distribution (SED) of the stars, we 
assume a Maraston (1998, 2005) single stellar population (SSP) 
parametrization. The age and metallicity of the SSP in the disk and 
bulge are fixed to 8 Gyr old and a solar metallicity {Z = 0.02), as 
obtained from diagnostic plots of several Lick indices (see Figures 
5 and 6 in Proctor et al. 2000). The age of the SSP is assumed to 
be the same in the stellar disk and bulge, since the boxy bulge is 
considered to be a pseudo bulge. The intensity of the stellar com- 
ponents is set by its V band luminosity (Lv). 

Earlier studies of the dust geometry in NGC4565 have sug- 
gested a ring-like structure for a large fraction of the dust (Engar- 
giola & Harper 1992; Ohta & Kodaira 1995; Neininger et al. 1996; 
Kormendy & Barentine 2010; Laine et al. 2010). Therefore, our 
dust model is composed of an exponential disk with a density dis- 
tribution as represented in Eq. 1 and a ring with a Gaussian distri- 
bution in radial direction and an exponential vertical profile 
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exp 
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where Ro, a and h z describe the radius of the ring, the radial 
dispersion and the vertical scale height. The content of the dust 
reservoir in NGC4565 is scaled by the dust mass in each of the 
dust components. The dust population is assumed to be uniform 
across the entire galaxy, consisting of a composition of dust parti- 
cles with a fixed grain size distribution. The abundances, extinction 
and emissivity of the dust mixture are taken from the Draine & Li 
(2007) model for the dust in our own Galaxy. 

With the different model components for the stars and dust 
in NGC 4565 and furthermore assuming a variable inclination for 
the galaxy, the degrees of freedom in the fitting procedure amount 
to 15 parameters. Our model does not account for the continuum 
emission from the AGN in the center of NGC 4565, since the low 
[Ne III] 15.55 /xm-to-[Ne II] 12.81 fim ratio (~ 1, Laine et al. 2010) 
suggest that the AGN continua are weak in this galaxy (Deo et al. 
2007). Table 2 gives an overview of the best fitting parameters ob- 
tained from the FitSKIRT algorithm. Since a sufficient number of 
photon packages were used to obtain this best fitting model (a pho- 
ton package is a collection of a number of photons, with the en- 
semble of all photon packages representing the radiation field in 
the galaxy), we assume an uncertainty of about 20% for all derived 
parameters (see also De Geyter et al. 2012). 

Exploring the compatibility of this model with optical data, 
we compare the observed g band image (see Figure 3, top row) 
and the major- and minor-axis g band profiles for several galacto- 
centric radii (see Figure 3, bottom row) with the modelling results. 
The resemblance between modelled and observed images, major- 
and minor-axis profiles suggests that our radiative transfer model is 
representative for the observed stellar emission and optical proper- 
ties of dust in NGC 4565. For the inclination, we find a best fitting 
value of rsj 87.25°, which is in agreement with earlier reported re- 
sults. Indeed, Hi observations of the gaseous disk in the past have 
constrained the inclination > 84° (Ruben 1991), with more recent 
Hi observations reporting an inclination angle of 87.5° (Zschaech- 



ner et al. in prep.). From a similar radiative transfer study, Alton et 
al. (2004) estimated an inclination of 88°. 

From the stellar distribution in this radiative transfer model, 
we infer a bulge-to-disk ratio of B-to-D ~ 0.47 in the V band. The 
model indicates a total dust mass of Ma ~ 1.0 x 10 8 M of which 
about 60 % resides in an exponential disk, with the remaining 40 % 
distributed in a broad ring structure. The scale size of the exponen- 
tial dust disk Qir ~ 223", h z ~ 2.5") agrees well with the results 
obtained from a similar radiative transfer modelling procedure in 
Alton et al. (2004) (h R ~ 200", h z ~ 4.3"). The diffuse dust com- 
ponent in Alton et al. (2004) has a face-on optical depth r v ~ 0.63, 
which corresponds to Md ~ 9.7 x 10 7 M©, when assuming the 
same extinction cross-section for the Draine & Li (2007) dust com- 
position applied in our dust model. 

Figure 4 represents the spectral energy distribution as mod- 
elled by our radiative transfer code SKIRT (solid, black line). Avail- 
able flux measurements were overlaid on this SED model to de- 
scribe the observed spectral energy distribution. An overview of 
flux measurements for NGC 4565 is given in Table 1. Less accu- 
rate flux measurements were considered non-relevant for this anal- 
ysis and were omitted from this plot (see Section 2). Although 
our model shows great resemblance with the g band image and 
major- and minor axis profiles, the spectral energy distribution 
in Figure 4 shows that our model fails to reproduce the emission 
from NGC 4565 in the UV and MIR/FIR wavelength domains. The 
higher observed UV and MIR emission from NGC 4565 indicates 
the presence of a young stellar population in this galaxy, which is 
not yet accounted for by our model. 

3.2.2 Model 2: standard model with star formation 

In an attempt to reproduce the ultraviolet and infrared emission ob- 
served from NGC 4565, we add a star formation component to the 
model constructed in Section 3.2.1. Rather than using an empiri- 
cal star formation template (e.g. MacLachlan et al. 201 1), we adopt 
starburst templates from the library of pan- spectral SED models 
for the emission from young star clusters with ages < 10 Myr 
presented in Groves et al. (2008). The SED templates in Groves 
et al. (2008) were generated from a one-dimensional dynamical 
evolution model of Hll regions around massive clusters of young 
stars, the stellar spectral synthesis code Starburst 99 (Leitherer et 
al. 1999) and the nebular modelling code MAPPINGS III (Groves 
2004). The Starburst99 templates used in these models correspond 
to an instantaneous burst with a Kroupa (2002) broken power-law 
IMF. The main parameters controlling the shape of the emission 
spectrum from these young stellar clusters are the metallicity (Z) 
of the gas, mean cluster mass (M c i), age and compactness (C) of 
the stellar clusters, the pressure of the surrounding ISM (Po) and 
the cloud covering fraction (/pdr). The metallicity is chosen similar 
to the metallicity of the old stellar population in NGC 4565 (Z ~ 
1 Z©). To eliminate the age parameter, we use age-averaged tem- 
plates obtained from averaging spectra over 21 cluster ages, rang- 
ing from 0.01 to 10 Myr in steps of 0.5 Myr. For the mean cluster 
mass and ISM pressure, we assume fixed values M c \ ~ 10 5 M© 
and Po/k ~ 10 6 cm -3 K. Those approximations are justified since 
more massive star clusters can be simply thought of as the superpo- 
sition of several individual clusters. A variation in the ISM pressure 
furthermore mainly affects the nebular emission lines rather than 
altering the shape of the emission spectrum (Groves et al. 2008). 
For the cloud covering fraction /pdr, we will assume a value of 
fpDR = 1 to describe the emission from heavily obscured star- 
forming regions in NGC 4565. The choices for M c \ and Polk leave 
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Figure 3. Model 1: the observed (left) and modelled (right) g band image for NGC4565 on the top row. The bottom row shows from left to right the major 
axis g band profiles at radial distances z=0, 10 and 30" from the center of NGC 4565, respectively. The rightmost panel represents the minor axis profiles at 
R=0. The observed profiles are indicated in black, with the modelled emission color-coded in green. 
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Figure 4. Model 1: the modelled SED as obtained with the SKIRT model (solid black line), consisting of merely an old stellar population and neglecting the 
contribution from young stars and/or a clumpy dust distribution, overlaid with the observed fluxes (see Table 1). 



us to explore the effect of the compactness parameter C on the 
shape of the SED for values of log C ranging between 4.5 and 6.5. 
With the compactness parameter controlling the position of the far- 
infrared emission bump by shifting it in peak wavelength without 
significantly altering the width of the emission curve, the parame- 
ter C is closely related to the dust temperature (Groves et al. 2008). 
Although the emission at 24 /im might still be dominated by hot 
dust emission predominantly heated by young stars of ages < 10 
Myr, the contribution of diffuse dust emission is non-negligible at 
far-infrared and sub-millimeter wavelengths. Therefore, it is diffi- 
cult to constrain the shape of the emission spectrum from hot dust 
grains and thus to determine the appropriate value of log C in the 
models. With values of 6.5 for log C describing extremely compact 
star-forming regions, which are characteristic for the high pressure 
ISM and/or massive young clusters in starburst galaxies, we assume 
a more moderate compactness factor of log C = 5.5. A value of log 
C = 5.5 is furthermore consistent with the dust temperature Td = 30 
K found for the warm dust component heated predominantly by 



stellar emission originating from star-forming regions (Engargiola 
& Harper 1992). 

In addition to these young stellar clusters, we have also in- 
cluded a component of ultra-compact Hll regions in our model. 
Those ultra-compact Hll regions are typically used to simulate the 
conditions during the earliest stages of the cluster lifetime (< 10 6 
yr) when the newly formed stars are likely still buried in their sep- 
arate birth clouds and a typical star cluster consists of an ensem- 
ble of individual ultra-compact Hll regions, rather than providing 
a single source of radiation originating from the smoothed emis- 
sion from all stars in the stellar cluster. The templates for these 
ultra-compact Hll regions were constructed in a similar way as for 
the dusty Hll regions with their surrounding PDRs, but each birth 
cloud is modelled individually and the template for the entire com- 
pact Hll regions is constructed by co-adding the individual SEDs 
(see also Groves et al. 2008 for a more thorough description of 
the models). The SED template for the ultra-compact Hll is jointed 
with the SED template for the Hll regions with fpDR = 1, under 
the assumption that they produce new stars at the same rate. This 
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joint SED template will be scaled until the observed emission at 
24 /im is reproduced by our model. 

Emission spectra for the young stellar population with ages 
between 10 and 100 Myr are taken from the Starburst99 library 
(Leitherer et al. 1999). The Starburst99 templates were constructed 
for an instantaneous burst with a solar metallicity, a burst mass 
M c i = 10 6 M and a Kroupa (2002) broken power-law IMF, which 
is consistent with the stellar emission spectra used to compute the 
SED of stellar clusters younger than 10 Myr in Groves et al. (2008). 
With the young stellar clusters < 10 Myr emitting predominantly in 
mid-infrared wavebands, the non-ionizing radiation of the evolved 
stellar population up to ages 50-100 Myr will dominate the UV 
emission spectrum in NGC4565 (Calzetti et al. 2005). 

Although most details on the exact location of star-forming 
complexes within the disk vanish along the line of sight, we try to 
recover the average geometry of unobscured and dust-enshrouded 
star-forming regions from the emission profiles along the major 
axis of NGC4565 in the GALEX NUV and MIPS 24 /im wave- 
bands (see Figure 5, top and bottom panel, respectively). By allo- 
cating a certain distribution to the star-forming regions, the proper 
amount of attenuation due to the foreground screen of diffuse dust 
is applied on every specific location within the galaxy. Indeed, due 
to the high inclination angle the different star-forming regions in 
NGC 4565 will experience wide ranges in obscuration depending 
on whether they are located on the near- or far end of the galaxy's 
disk with respect to the line-of- sight. From the major axis NUV 
emission profile, we derive a ring-like geometry for the unobscured 
star formation component, modelled as a ring (see Equation 4) with 
central radius R ~ 290" and width a ~ 72". The scale height 
of the ring of young stars is assumed to be identical to the scale 
height of the stellar disk populated with the old stellar population 
(h z ~ 8.2"). From the MIPS 24 /im image, it becomes immedi- 
ately evident that a component of obscured star formation resides 
in a narrow ring coinciding with the CO molecular ring centered 
at 80-100" (Sofue & Nakai 1994; Neininger et al. 1996). The pa- 
rameters for the ring (see Equation 4) of embedded star formation 
are chosen similar to the molecular CO ring with central radius 
at R ~ 90", width a ~ 3.4" and scale height h z ~ 1.5". Aside 
from this concentration of embedded localized sources distributed 
in a ring, a more extended component of obscured star formation 
seems also present in NGC 4565. To reproduce the more extended 
24 /im emission in the major axis profile, obscured star-forming 
complexes are distributed in the ring harboring the unobscured star 
formation component and the dust ring containing ~ 40% of the 
diffuse dust component in NGC 4565. From the 24 /im major axis 
profile (see Figure 5, bottom panel), the presence of a central dust 
disk becomes evident. This central dust concentration is modelled 
as an exponential disk with radial and vertical scale lengths of 5" 
and 1", respectively. The distribution of this central dust reservoir 
(M d ~ 3 x 10 4 M ) is truncated at a radius of 20". 

When fitting the UV and 24 /im emission of NGC 4565, we 
require an average star formation rate in the outer regions of the 
disk (with most of the young stars between galacto-centric radii of 
200" and 400") of 3.1 M yr" 1 for the 10 to 100 Myr old stel- 
lar population. The heavily obscured star-forming regions of ages 
< 10 Myr seem to reproduce the observed 24 /mi emission when 
assuming an average star formation rate of 2.2 M yr -1 over the 
past 10 Myr. The dust-enshrouded star-forming clumps in the nar- 
row ring produce young stars at a rate of SFR ^0.6 M yr -1 , with 
a somewhat higher productivity for the dense clouds distributed in 
the ring of unobscured star formation (SFR ^0.8 M yr -1 ) and 
the diffuse dust ring (SFR ^0.8 M yr -1 ). This obscured star for- 
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Figure 5. Model 2: NUV (top) and MIPS 24 /im (bottom) major axis pro- 
files. The observed profiles are indicated in black, with the modelled emis- 
sion color-coded in red. 



mation rate inferred from the models is perfectly consistent with 
the SFR - 2.0 M yr" 1 calculated in Wu & Cao (2006) from the 
far-infrared luminosity Lfir in NGC 4565, corrected for the contri- 
bution from an old stellar population (Satyapal et al. 2005). Also 
monochromatic SFR tracers such as the 24 /im emission result in 
similar estimates for the SFR ~ 2.0 M yr -1 (Rieke et al. 2009). 
We argue furthermore that the star formation rate of stars between 
10 and 100 Myr in our model (SFR - 3.1 M yr" 1 ) is in close 
agreement with the SFR ~ 2.7 M yr -1 calculated from the FUV 
luminosity Lfuv for NGC 4565 and relying on the SFR calibration 
in Salim et al. (2007). 

Accounting for all young stellar objects in our model, we de- 
rive an average star formation rate of 5.3 M yr -1 over the past 
100 Myr. The small decrease in the star formation activity during 
the recent 10 Myr is in agreement with a delayed exponential law 
describing the star formation history in other late-type spiral galax- 
ies (Boselli et al. 2001; Gavazzi et al. 2002). The total star forma- 
tion activity in our model (SFR ^5.3 M yr -1 ) is furthermore 
consistent with the total star formation rates derived from a com- 
bination of unobscured and obscured star formation tracers (e.g. 
SFRpuv+24/xm ~ 3.2 M yr" 1 (Zhu et al. 2008), SFR FU v + tir ~ 3.9 
M yr" 1 (Buatet al. 2005)). 

This model with both unobscured and obscured star forma- 
tion is able to account for the young stellar emission at UV/mid-IR 
wavelengths (see images on the first and third row of Figure 6, re- 
spectively) and remains in agreement with the optical constraints 
(see images on the second row and plots on the last row of Fig- 
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ure 6) upon reducing the luminosity of the old stellar population 
in the disk by a small factor to correct for the emission from those 
young stars in optical wavebands (Lv = 4.0 x 10 10 L©,v). The to- 
tal dust mass associated with the young stellar clusters (< 10 Myr) 
in our models is Md 3.6 x 10 7 Mq, which increases the total dust 
mass in NGC4565 to M d 1.4 x 10 8 M . The geometry for the 
embedded, localized sources in NGC4565 is able to account for 
the 24 fim emission in NGC 4565 as observed from the major axis 
profile (see Figure 5, bottom panel) or the MIPS 24 /im image (see 
images on the fourth row of Figure 6). The geometry for the less 
obscured young stars with ages > 10 Myr is capable of reproducing 
the emission profile along the major axis (see Figure 5, top panel) 
as well as the emission originating from the entire galaxy at NUV 
wavelengths (see Figure 6, third row). The peak at R ~ 420" in the 
NUV major axis profile belongs to the residual star formation in the 
warp structure on the northwestern side of the galaxy's disk and is 
not taken into account in our model. The stellar clusters do have a 
more clumpy distribution throughout the disk than was accounted 
for by our model. 

More critical is the inconsistency in the FIR/submm, where 
our radiative transfer model continues to underestimate the ob- 
served emission spectrum for NGC 4565 (see Figure 7). Even 
though the model accounts for the dust heating provided by star 
formation, the modelled dust emission underestimates the observed 
far-infrared and submillimeter emission by a factor of ~ 3-4, in cor- 
respondence to dust energy balance studies for other nearby edge- 
on spirals. 

At mid-IR wavelengths, our model also underestimates the 
observed emission originating mainly from PAHs in those wave- 
bands by a factor of ^ 2 (see Figure 7). This dissimilarity might 
on the one hand be a consequence of the fixed PAH template used 
in the models to construct the emission spectra reported in Groves 
et al. (2008). Based on the average variation by a factor of 2 ob- 
tained from a study of the mid-IR emission for a large sample of 
nearby galaxies (Smith et al. 2007), we should interpret this off-set 
as a diagnostic for the specific ISM conditions in NGC 4565 rather 
than attributing it to a deficiency of star formation in our model. 
With the PAH template in Groves et al. (2008) being optimized 
to fit the Spitzer IRS observations of two interacting galaxy pairs 
(NGC 4676 and NGC 7252), small deviations in ISM conditions 
and interstellar radiation field hardness might explain the difference 
in observed PAH emission for NGC 4565. On the other hand, the 
diffuse emission from cold dust at FIR/submm wavelengths has not 
yet properly been accounted for by our radiative transfer model and 
might be responsible for an important fraction of the PAH emission 
in galaxies (e.g. Haas et al. 2002; Bendo et al. 2008). 



4 DISCUSSION 

4.1 Dust energy balance 

Although in agreement with the observed properties of stars and 
dust at optical wavelengths and the young stellar emission in UV 
and MIR wavebands, the model from section 3.2.2 underestimates 
the observed dust emission in the far-infrared and (sub)millimeter 
wavelength regime. Whereas the dust emissivity of dust grains 
could influence the shape of the SED at some wavelengths (Al- 
ton et al. 2004), the inconsistency in the integrated dust energy 
balance can not be solely caused by a grain population with dust 
emissivities diverging from Galactic values at longer wavelengths. 
Indeed, the multi-frequency sampling of the dust SED allows us de- 
rive a true discrepancy in the integrated dust energy balance, which 



Table 2. Best fitting model parameters for the stellar (exponential 
disk + outer bulge) and dust components (exponential disk + ring) in 
NGC 4565 as obtained from FitSKIRT. 



Geometry 




87.25 


Stellar disk 


l r//n 

[ ] 

h z ["] 

Ly [L^v] 1 


102 
8.2 
4.9 x 10 10 


Outer bulge 


n 

r ef f ["] 

q 

Ly [L©,v] 


2.12 
20.2 
0.744 
2.0 x 10 10 


Dust disk 


h R ["] 
h z ["] 
M d [M ] 


223 
2.5 
6.3 x 10 7 


Dust ring 


R["] 
h z ["] 
*["] 

M d [M ] 


146 
4.3 
41.2 
4.0 x 10 7 



lr The luminosity of the stellar disk is later adjusted to Ly = 4.0 x 10 10 
Lq^v to correct for the emission from young stars in the V band when 
supplementing our model with star formation. 

is impossible to solve by assuming divergent FIR emissivities due 
to a different grain composition and structure compared to the typi- 
cal Galactic dust reservoir. Although we can not exclude a possible 
variation in dust emissivity at submm/mm wavelengths (diverging 
from the assumed emissivity in the Draine & Li 2007 dust model), 
we argue that this discrepancy is likely due to a clumpy morphol- 
ogy of the dust in NGC 4565 unresolved in the currently available 
dataset of far-infrared/submillimeter observations for NGC 4565. 
Due to the high density in those clumps, they hardly contribute 
to the attenuation of stellar light and can therefore not be distin- 
guished from optical observations. 

Due to the edge-on inclination of NGC 4565 and the current 
limitations in resolving power of infrared instrumentation, we are 
not capable of constraining the exact location and size of those dust 
clumps and rather have to remain content with information on the 
average mass and temperature of this clumpy dust reservoir. Those 
dust clumps either correspond to quiescent clouds where star for- 
mation has not (yet) been initialized. On the other hand, this extra 
dust reservoir might reside in the outer shells of star-forming com- 
plexes where the heating by the young embedded objects becomes 
negligible due to the high optical depth of the dust shells in the 
immediate vicinity of the stellar source. 

To estimate the amount of dust residing in clumps, we fit the 
DustEm dust model with a Draine & Li (2007) grain composi- 
tion to the residual fluxes after subtracting the dust emission in 
our model from the observed flux densities (IRAC 5.8, 8.0 /xm, 
WISE 12.1 /mi, IRAS 12 and 100 /mi, Kuiper 100 and 160 /mi, 
MIPS 160 /mi, SPIRE 250, 350, 500 /mi, IRAM 1.2 mm). Figure 
7 shows the resulting SED (black, solid curve) when supplement- 
ing our model with a dust reservoir of ^ 1.5 x 10 8 M at a tem- 
perature of Td ~ 16.9 K (green, dotted curve). The average tem- 
perature Td ~ 16.9 K obtained for the residual dust reservoir is 
consistent with the temperature of the diffuse dust in the disk of 
NGC 4565 {Td ~ 14-18 K) and the average colour temperatures 
derived in other nearby spiral galaxies at wavelengths longwards 
of 160 /xm (e.g. M81, M83 and NGC 2403, Bendo et al. 2012). 
At mid-IR wavelengths, this additional dust reservoir is capable of 
explaining the PAH emission in NGC 4565. Also at FIR/submm 
wavelengths, a massive cold dust component seems adequate to 
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Figure 6. Model 2 with on the top row: the observed (left) and modelled (right) NUV band (top row), g band (second row), 24 /im (third row) and 250 /im 
image (fourth row). The bottom row represents from left to right the major axis g band profiles at radial distances z=0, 10 and 30" from the center of 
NGC4565, respectively. The rightmost panel represents the minor axis profiles at R=0. The observed profiles are indicated in black, with the modelled 
emission color-coded in green. 




Figure 7. Model 2: the modelled SED as obtained with the SKIRT model (dashed purple line) overlaid with the observed fluxes (see Table 1). Model 2 
corresponds to model 1 supplemented with young star clusters (< 100 Myr). To estimate the dust mass residing in quiescent dust clumps in NGC4565, 
we fit the DustEm dust model with a Draine & Li (2007) grain composition (green, dotted line) to the residual fluxes for NGC 4565 in the far-infrared and 
submillimeter wavebands. The model with star formation and an additional dust reservoir distributed in clumps can explain the multi- wavelength observations 
for NGC 4565 (see black, solid curve). 
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explain the emission in those wavebands. At 60 and 70 /xm, our ra- 
diative transfer model however overestimates the observed fluxes. 
This inconsistency could be an indication for the star-forming com- 
plexes in NGC 4565 to be more compact than currently accounted 
for by our model. Due to the high inclination angle of NGC 4565, 
a direct identification of heavily obscured star-forming complexes 
within the disk of the galaxy and their individual size distribution 
and compactness is impeded. 

With these additional cold dust clumps, the total dust mass in 
NGC 4565 amounts to M d ~ 2.9 x 10 8 M of which one-third 
has a diffuse distribution throughout the disk and the remaining 
two-thirds resides in clumps. From those dense dust clouds, the 
majority (~ 80%) is not heated by UV-emitting sources while the 
remainder (~ 20%) host embedded star-forming complexes. The 
total dust mass Md ~ 2.9 x 10 8 M in our model is consistent with 
the dust reservoir accounted for by BLAST observations (Wiebe et 
al. 2009), but ^10 times higher than the dust content reported in 
Engargiola & Harper (1992) (2.7 x 10 7 M ) based on IRAS ob- 
servations up to 100 /xm. This confirms the need for longer wave- 
length observations to trace the entire dust reservoir in a galaxy 
(e.g. Bendo et al. 2010; Gordon et al. 2010; Bendo et al. 2012). 

The solution to the dust energy budget problem for NGC 4565 
accounting for a sizable fraction of the dust (~ 65%) distributed in 
compact clumps is in agreement with other energy balance studies 
of edge-on objects (e.g. Bianchi 2008; Baes et al. 2010; Popescu et 
al. 201 1; De Looze et al. 2012; Holwerda et al. 2012, etc.). 

4.2 Gas-to-dust ratio 

From the Hi mass M H i ^ 1.8 x 10 10 M (Zschaechner et al. in 
prep.) and the CO observations (M H2 ~ 2.9 x 10 9 M ) presented 
in Neininger et al. (1996) (both quantities are already converted to 
our adopted distance of D = 16.9 Mpc), we derive a total gas mass 
of ~ 2.9 x 10 10 M for NGC 4565 when applying a correction 
factor of ~ 1.4 to correct for heavier elements. Based on the total 
dust mass in our radiative transfer model (Md ~ 2.9 x 10 8 M ), 
we estimate an average gas-to-dust ratio of ~ 100 in NGC 4565, 
which is in fair agreement with the gas-to-dust fraction found in 
other late-type spirals (e.g. ~ 120 in our own Galaxy, Zubko et al. 
2004) and the dust-to-Hl mass fraction derived for Sb galaxies with 
normal Hi components (as opposed to Hi deficient objects) in the 
HRS sample of galaxies (Cortese et al. 2012a). 

4.3 Dust heating sources 

Investigating the dust heating in NGC 4565, we find that about 70% 
of the far-infrared luminosity is attributed to the absorption of stel- 
lar light from the old stellar population, with the heating for the re- 
maining 30% provided by the young stellar population. In this cal- 
culation, we assume that the high-density dust clumps distributed 
throughout the diffuse ISM in NGC 4565 are merely heated by the 
old stellar population, since the temperature of the dust (Td ~ 16.9 
K) in those dense dust cores does not suggest any link to localized 
embedded sources and is more closely related to the average tem- 
perature of the diffuse dust in the disk of the galaxy (Td ~ 14-18 
K). Those heating fractions are in agreement with similar studies 
of dust heating mechanisms in other nearby galaxies, where the old 
stellar population also accounts for an important fraction of the far- 
infrared radiation (e.g. Bianchi 2008; Bendo et al. 2010; Boquien 
et al. 201 1; Popescu et al. 201 1; Bendo et al. 2012). Similar results 
were obtained from accurate modelling of the absorption of stel- 
lar light as well as careful investigations of UV observations (Buat 



& Xu 1996; Chariot & Fall 2000; Witt & Gordon 2000; Cortese 
et al. 2008). Engargiola & Harper (1992) report a dust heating ef- 
ficiency of 54 and 46% for the young and old stellar population 
in NGC 4565, respectively. Considering that their IRAS observa- 
tions only covered a wavelength range from 12 to 100 /mi, it is not 
surprising that they underestimated the contribution from the old 
stellar population to the dust heating with the more evolved stars 
contributing the bulk of dust heating at longer wavelengths. 

Although the fraction of dust heated by the young stellar pop- 
ulation in NGC 4565 is overall small (~ 30 %), the contribution of 
hot dust emission from star-forming complexes is about 90% in the 
24 fim waveband, after which it decreases quickly to about 50% at 
70 fim, 12% at 160/xm, less than 10% in all SPIRE wavebands and 
below 5% at 1.2 mm. Those values are consistent with the decreas- 
ing contribution of star-forming regions for increasing wavelength 
in M81, M83 and NGC 2403 (Bendo et al. 2012), with a relative 
contribution of the young stellar population to the total dust heat- 
ing of more than 50% shortwards of 250 /im and dropping below 
30% in the SPIRE wavebands. 

This self-consistent analysis of the dust heating mechanisms 
in NGC 4565 based on a realistic radiative transfer model con- 
firms earlier studies reporting on a non-negligible fraction of the 
dust heating in spiral galaxies powered by the more evolved stel- 
lar population. With the young stellar population being responsi- 
ble for less than half of the total infrared dust emission originating 
from NGC 4565, we believe caution is needed when using the to- 
tal infrared emission in galaxies to trace the star-forming activity. 
Although a tight correlation was found between the star formation 
and the total infrared dust emission for large samples of nearby spi- 
ral galaxies (Devereux & Young 1990; Devereux et al. 1995; Buat 
& Xu 1996; Kennicutt et al. 2009), this relation might be the result 
of an indirect link between the star-forming activity in galaxies and 
the total-infrared emission in galaxies, tracing the surface density 
of gas in the interstellar medium. In the latter case, the SFR-Lfir 
correlation is governed by the Schmidt law (Schmidt 1959; Ken- 
nicutt 1998), relating the surface density of gas and star formation 
in galaxies, rather than the dust heating provided by star-forming 
regions. 



5 CONCLUSIONS 

We present a full radiative transfer analysis of the edge-on spiral 
galaxy NGC 4565, accounting for the absorption and scattering of 
stellar light by dust grains and its thermal re-emission at infrared 
wavelengths. From a radiative transfer model fitting procedure to 
the optical SDSS g band image, we determine the best fitting pa- 
rameters for the old stellar population and diffuse dust component 
in NGC 4565. To account for the observed UV and mid-infrared 
emission in NGC 4565, we supplement our model with a young 
stellar population of age < 10 Myr (SFR ~ 2.2 M yr _1 ) and rang- 
ing from 10 to 100 Myr (SFR ~ 3.1 M© yr -1 ). The distribution of 
star-forming complexes within the disk of the galaxy is constrained 
from major axis NUV and 24 fim emission profiles. 

Even though this young stellar population provides an addi- 
tional power source for dust heating and hereby boost the emission 
at mid-infrared wavelengths, the emission observed at wavelengths 
longwards of 100 /mi remains underestimates by a factor of 3-4. 
This inconsistency in the dust energy budget of NGC 4565 sug- 
gests the presence of a sizable fraction (two-thirds) of the total dust 
reservoir (Md ~ 2.9 x 10 8 M©) in a clumpy distribution with no 
associated young stellar sources. The distribution of those dense 
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dust clouds would be such that they remain unresolved in current 
far-infrared/submillimeter observations and hardly contribute to the 
attenuation by dust at optical wavelengths. The contribution from 
a grain composition with dust emissivity properties diverging from 
the typical Galactic standards can however not be discarded as a 
factor of influence on the high far-infrared/submillimeter fluxes. 
The majority of dust heating in NGC4565 is provided by the old 
stellar population (70%), with the remaining dust heating powered 
by localized embedded sources in NGC 4565. 

The results from this detailed dust energy balance study in 
NGC 4565 accords with similar analyses of other edge-on spirals, 
concluding that a significant fraction of the dust is distributed in 
such a way that the influence on the overall extinction of a galaxy 
is negligible but does emit strongly in far-infrared/submm wave- 
bands. 
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